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Abstract 

We study the evolution of regions of space with various initial field values for a simple 
theory that can support hybrid inflation. Only very narrow domains within the range of 
initial field values below the Planck scale lead to the onset of inflation. This implies a severe 
fine tuning for the initial configuration that will produce inflation. 
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One of the most attractive feature of hybrid inflation |1| is that it can take place for values 
of the inflaton field below the Planck scale This permits the construction of realistic field 
theoretical models without the need of a description of quantum gravity. Several such models 
have been proposed within the context of supersymmetry Inflation requires the presence 

of an almost flat direction in field space with non-zero vacuum energy. A very small slope along 
this direction results in the slow rolling of the inflaton field (j). For hybrid inflation, the slow 
rolling occurs along a valley of the effective potential. It ends when the valley turns into a ridge 
and the slow-roll regime terminates, owing to the growth of fluctuations of fields orthogonal to 
the inflaton field. Typically this point of instability corresponds to a value (pins of the inflaton 
below the Planck scale. 

The question we would like to address is whether generic initial values of the fields below 
the Planck scale can lead to the onset of inflation. Clearly, there is an area where this is likely, 
the almost flat direction. Our aim is to explore the whole field space below the Planck scale and 
establish whether inflation is a frequent occurrence or a rare event. Similar investigations have 
already been performed in refs. [10, 0]. Our main observation is that the field range that can 
lead to inflation is much more constrained than what one would naively infer from refs. |1C, p]] , 
unless a significant fine tuning of the fields is imposed within the space regions that eventually 
inflate. 

The onset of inflation requires a region of space, the size of a few Hubble lengths, where the 
fields take almost constant values, so that the gradient energy density is negligible compared 
to the potential energy density. The earliest time at which one could start talking about such 
regions of space is when the Planck era (during which quantum gravitational fluctuations dom- 
inate) ends and classical general relativity starts becoming applicable. The initial field values 
within each region are expected to be of order Mpi, or larger if the couplings of the theory are 
small, as in chaotic inflation jl^. In this work we concentrate on field values below Mpi for 
two reasons: a) Hybrid inflation does not require extremely small couplings, b) In general, the 
field-theoretical models for which a potential can be reliably computed assume an expansion in 
powers of cp/Mpi. 

Inflation could start at the end of the Planck era if the fields take values in the parts of 
the potential that are able to support it. However, there is also the possibility that the fields 
will evolve, from some initial values that do not give inflation, to different values that do. We 
consider all possible field values, assuming an initial degree of homogeneity over a few Hubble 
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lengths, and follow the evolution by integrating the classical equations of motion. In this way 
we determine the range of field values that eventually lead to inflation. We show that, for many 
realistic models, the initial fields must be extremely homogeneous if their initial values are below 
the Planck scale. For example, for the prototype model of hybrid inflation (see below), one of the 
fields must be homogeneous with at least ~ 10~^-10~^ accuracy. It is very difficult to calculate a 
probability distribution for the initial field configurations, and thus have a quantitative estimate 
of the probability of inflation to occur. This requires a reliable description of Planck scale 
dynamics. For this reason we do not address the question of whether such a homogeneous initial 
state is probable or not. Our results have a much looser interpretation, such as indicating how 
"natural" hybrid inflation is. 

We consider the simplest potential that can support hybrid inflation Q 



with positive mass terms and couplings. It possesses the symmetries cj) — > —(j) and a — > —a. 
A renormalizable term A'(/>^ is permitted by these symmetries and should naturally appear. 
However, neglecting such a term can be justifled if the hybrid inflationary scenario is embedded 
in the context of supersymmetry and an appropriate superpotential is chosen. The term m^(/>^/2 
is also absent in supersymmetric models, and the small slope along the inflationary trajectory 
is generated by non-renormalizable (such as logarithmic) terms that are present when quantum 
corrections are taken into account ^, P]. 

We use the simplest choice of eq. (|^) for our study. Our results do not depend on whether the 
slope along the inflationary trajectory is generated by the term jl or radiative corrections. 
In the concluding paragraphs we discuss a simple superpotential to which our results apply 
directly. We point out, however, that our analysis must be repeated for the determination of 
the appropriate initial conditions in supersymmetric inflationary models with more complicated 
superpotentials. 

The potential of eq. (|lj) has a valley of minima along the ^ axis [a = 0) for flxed 0^ > ^f^^. 



There is a series of maxima along the </> axis (a = 0) for flxed (p'^ < cff^s- The absolute minima 
of the potential are located at ((/> = 0, cr = zizM). We assume that a Robertson- Walker metric is 
a good approximation for the regions of space with uniform flelds that we are considering. The 





where 
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evolution of the fields is given by the standard equations 
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with Mpi = 1.2 X 10^9 GeV. 

Inflation occurs along the (p axis, where the potential is 

V{4>) = ^AM^ + im2</.2. 



(6) 



In the limit where the second term in the r.h.s. of the above equation dominates, we obtain 
the prototype model of chaotic inflation. This takes place for values of the field cf^ S> Mp^. 
The COBE observation of the cosmic microwave background anisotropy leads to the constraint 
m ~ 10^^ Mpi |12|. As we have already explained, we are interested in values of the field 



02 <, Mpi, for which the first term in the r.h.s. of eq. 
occurs for 

The "slow-roll" parameters |13| | are given by |Q 
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In the interval 
for the parameter range of eq. (0). 

The number of e-foldings that occur between the moment that cosmologically interesting 
scales (of the order of the present Hubble distance) leave the horizon and the end of infiation is 
~ 60. We denote by <j)QQ the value of the inflaton fleld at horizon crossing for these scales. In 
our model, inflation terminates when rolls beyond the instability point (pins- In the parameter 
range of eq. (^), the relat ion between (pQo a^nd. (pins is 0] 
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with r] given by eq. (^). The inflationary prediction for the spectrum of adiabatic density 

perturbations at cosmological scales |jl^, ^ is 

2 ^ 32 Feo 1 ^, Stt X^M^^ 
" 75 M|,; eeo ^ 75 Ml^mVio ' 

Comparison with the value 6h = 1-94 x 10~^, deduced from the COBE observation of the cosmic 

microwave background anisotropy, leads to the constraint p| 

M „.„ io„ _i/in,_i/c;/ m \^/^ 
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For couplings that are not many orders of magnitude smaller than 1 and a sufficiently small 
"slow-roll" parameter rj, we obtain a simple relation between the scales M and m 

M f m \'^/^ 

-—-0.13 -— . 13 

Mpi \MpiJ 

The above relation holds as long as our assumption that the first term in the r.h.s. of eq. (P) 
dominates is valid. As we have already mentioned, when the second term dominates, this model 
reduces to the prototype model of chaotic inflation, for which the COBE observations imply 
m = 10~^ Mpi. We expect, therefore, that eq. (^) holds as long as m ^ 10-^ Mpi. This leads 
to the constraint M ^ 5 x 10~'*Mp;. A detailed analysis Q verifies this conclusion and predicts 
a more accurate bound Af ^ 2 x lO^^Mp;. A similar bound applies to all models of hybrid 
inflation with couplings of order 1. On general grounds one expects the inflationary energy scale 
V^^^ to be at least two or three orders of magnitude smaller than the Planck scale [p^]: 

In order to determine the range of field values that lead to inflation, we start by considering 
a region of space, the size of a few Hubble lengths, in which the fields are homogeneous. We 
determine the evolution of this region by integrating eqs. (|^~(^ numerically. We consider field 
values < 0, cr < Mp;, and couplings X = g = 1. The choice of couplings has a small effect 
on our conclusions, as long as they are not much smaller than 1. Since we are interested in 
the necessary conditions for the onset of inflation, we neglect the term —m?(j) in sq. (^ during 
the numerical integration. We assume that inflation occurs if the evolution stops on the <j) axis 
beyond the instability point. If there is a significant slope along the (j) axis, the range of field 
values that lead to inflation is smaller than our predictions. 

In fig. 1 we display an array of initial field values in Planck units for a theory with M = 
IQ-^Mpi. The dark hexa gons denote initial conditions that lead to inflation along the positive 
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4> axis. The trajectories are short, with the (j) field moving directly to its final value, while the 
a field performs damped oscillations around zero. The white hexagons denote non-inflationary 
trajectories that end up at the minimum at (0 = 0,0" = M). The white triangles denote 
trajectories that end up at the minimum at = 0,cj = —M). The white squares at the 
beginning of the (p axis denote non-inflationary trajectories that start below the instability point 
{(j) = \J\/gM, (7 = 0) and end up at one of the two minima through the growth of fluctuations 
of the (j) field. We observe a large continuous range of initial conditions that lead to the onset 
of inflation. This range starts at the instability point. 

In fig. 2 we display a similar array of initial field values for M = 3 x IQ-^Mp/. The black 
triangles indicate initial values that lead to inflation along the negative (j) axis. We observe that 
the large continuous range of inflationary initial values has moved towards the Planck scale. A 
thin strip remains along the (p axis that ends at the instability point. Trajectories starting in the 
rest of field space rarely lead to the onset of inflation. Both fields oscillate around zero before 
the system relaxes either at one of the minima or in the valleys along the (j) axis. In fig. 3 we 
present a magnification of fig. 2 around the origin, which demonstrates the size of the various 
areas that lead to inflation. The largest of these has a size ~ 2-3 x lO^^Mp;, but the typical 
size is ~ l{)~^Mpi. The thin strip around the 4> axis has a width ~ 5 x Another 
magnification of fig. 2, away from the origin, is presented in fig. 4. Inflationary strips of width 
~ 10~^ Mp; are visible. 

Another array of initial field values, for M = 10~^Mp;, is presented in fig. 5. The situation 
is similar to that of fig. 2, but the large inflationary range has moved further towards the Planck 
scale. Figures 6 and 7 are magnifications of fig. 5 that demonstrate the size of the various areas. 
The inflationary strips of fig. 6 have a width ~ 1-2 x lO^^Afp^. Both fields oscillate many times 
before the evolution stops. The strip along the cf) axis has a width ~ 4 x 10~^Mp;. 

The shrinking of the areas of field space that lead to inflation is due to the decrease of H 
during the last stages of the evolution as M is lowered. This reduces the magnitude of the 
friction term in eqs. (§) and (^. Consequently, it takes longer for the system to settle down 
either at the minima or in the valleys along the (j) axis. Many oscillations of the fields are 
required before most of the initial potential energy is dissipated through expansion. This results 
is increasing sensitivity to the initial conditions. A slight change of the initial field values causes 
the system to move off an inflationary trajectory and end up at one of the minima instead. For 
M <, IQ-'^Mpi the largest area that leads to inflation, for initial field values not very close to 
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the Planck scale, is the strip along the (f) axis. There, the evolution is rather simple, with only 
the a field oscillating and cj) moving directly to its final value. 

We can obtain an estimate of the width of this area if we consider eq. (Q) with = 
and substitute for cr^ its average value u^ms — o'^/2. There are two time scales characterizing 
the solutions of this equation. The first one is related to the friction term and is given by 
t"! = 3H/2. For cr < M and (p/M < ./X/g M/a, we have 



The other time scale is obtained if we neglect the friction term and consider the oscillations of 
the (p field. One-fourth of the period gives the typical time for the system to roll to the origin 
and away from an inflationary solution. It is given by 

vr 1 

tosc = (16) 



/2ga 

Inflation sets in if tosc > tn, which gives 



For comparable couplings A, g, our assumptions for the derivation of the above bound break 
down when either M ~ Mpi or ^ y/AjZi^Mpi ~ 0.2Mpi. This approximate analysis is verified 
by the results in figs. 3 and 7. We have also checked eq. ( p^ numerically for M = \{)~^Mpi. For 
this value, the large continuous range of inflationary initial conditions moves to (j) field values 
above Mpi. 

The implications for the initial configuration that will lead to the onset of inflation are severe. 
For initial field values below Mpi and M ^ 10^'^ Mp;, the most favourable area of inflationary 
initial conditions is around the (j) axis. Throughout a region of space the size of the order of the 
Hubble length (which initially is ~ Mp^), the spatial variation of the a field in Planck units 
cannot exceed the r.h.s. of eq. (p!7|). If this condition is not satisfied, the fields in different 
parts of the original homogeneous space region will evolve towards very different values. In one 
part they may end up in the valleys along the (p axis, while in another they may settle at the 
minima of the potential. Before inflation, the size of space regions shrinks compared to the 
Hubble distance (for a quantitative discussion, see below). As a result, large inhomogeneities 
are expected at scales smaller than ~ H^^ when the evolution of the fields finally slows down. 
These will prevent the onset of inflation. We conclude that, for realistic inflationary models 
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(that require M ^ lO^^Mp;), a significant fine tuning is required for tlie initial configuration 
that will lead to inflation. 

The argument that we have presented is not rigorous. In principle, one should integrate the 
field equations for a configuration with a small spatial dependence and observe the growth of the 
perturbation However, we believe that the approximation of considering the initial space 
region as being composed of a few smaller regions with slightly different constant field values 
captures the essence of the problem. For this reason, we expect that a more rigorous treatment 
will not invalidate our qualitative conclusions. 

Up till now we have not considered the kinetic energy density of the fields. We have assumed 
that the fields evolve from a homogeneous value with zero time derivatives. This is not a good 
approximation when the potential is much smaller than Afp^. A more reasonable approximation 
is to fix the total initial energy density at Mpj, and consider a kinetic term equal to the difference 
between the total energy density and the potential. However, this does not fix the direction of 
the initial time derivative of the fields. If this direction is uniform throughout the initial space 
region, we do not expect a qualitative change in the behaviour we have observed previously 
In fig. 8 we display an array of initial field values in Planck units for a theory with M = 10~^ Mp/. 
The initial energy density has been fixed sA, E = 10~^Mpj. The initial time derivatives of the 
fields have been chosen equal and positive. We observe a strong similarity between figs. 6 and 
8. This is explained by the fact that the initial time derivative quicly causes the displacement 
of the fields from the </> axis into an area where 0, a ~ O.lMp;. There, the pattern of fig. 6 is 
expected. The larger inflationary strip of flg. 7 has disappeared. If the initial energy density is 
taken E ~ Mp^, the fields are displaced to regions (f),a ^ Mpi, where the inflationary strips are 
narrower than the ones depicted in fig. 8. We conclude that the presence of an initial kinetic 
term imposes much stronger constraints on the initial conditions than that of eq. (0)- There 
is the exception of an initial time derivative along the (p axis. However, this possibility leads to 
field values larger than Mpi. The variation of the direction of the time derivative within the 
initial space region results in more irregular patterns than that of fig. 8, and strengthens the 
constraints further. 

We have also followed the evolution of the scale factor R relative to the Hubble parameter, 
starting from an initial value Rq ~ Hq^. For the initial conditions of fig. 8 that lead to infiation, 
we find that at the onset of inflation R^^ is typically smaller than H^^ by a factor of order 10. 
For initial energy densities E ~ Mpi, this factor becomes of order 100. This implies that the 
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initial homogeneous region should extend far beyond a few initial Hubble lengths for this region 
to inflate. This factor becomes of order 1 if the initial time derivative of the fields is small and 
the fields start near the 4> axis and move quickly on it. 

We have based our discussion on the integration of the classical equations of motion. The 
presence of quantum (or thermal) fluctuations of the fields around their average values can be 
treated as a stochastic noise. We expect that this contributes an additional destabilizing factor, 
as it can easily knock the fields off an inflationary trajectory and enhance the irregularity of 
patterns in figs. 1-8. 

We conclude that hybrid inflation with the potential of eq. (||) is not "natural" . An extremely 
homogeneous initial field configuration must be assumed for initial field values below Mpi. For 
realistic models with M ^ 10~^Mp;, the a field must be homogeneous with at least ~ 10~^-10~^ 
accuracy. This bound may be strengthened by several orders of magnitude if the initial fields 
have a significant time derivative. These constraints do not apply if the initial values of the fields 
are larger than Mpi. One could imagine that inflation starts in this region and continues as 
the system rolls down the valley on the (j) axis. It is not possible, however, to construct reliable 
field-theoretical models in the field range above Mp;. 

Our results have direct application to the simplest supersymmetric model that supports 
hybrid inflation. It is described by the superpotential [0, 0] 



W = ^ 



Here ^ and ^ are chiral superfields, for which we assume canonical kinetic terms. The above 
superpotential is the only renormalizable one consistent with a continuous U{1) /^-symmetry 
under which W — > e^^W, $ — > e*^<I>, — > The superfields ^' and ^ transform under 

an internal gauge symmetry, while the $ superfield is a gauge singlet. The gauge symmetry is 
spontaneously broken at the scale [il\f\. The scale of inflation is set by // and is constrained 
by the COBE observations to be two or three orders of magnitude below Mpi. 

The potential of the scalar components of the superfields (for which we use the same notation) 
is given by 



V 



/i^ + A^'^' +A^|$r l^'r + l^'l +Z)-terms. (19) 



Vanishing of the D-terms is achieved along the D-flat directions where |^'| = |^'|. If the fields 
are assumed real, the potential along the D-flat directions is completely analogous to that of eq. 
(11]) but without a mass term for <I>. The slope along the inflationary trajectory can be generated 



either through the addition of a soft supersymmetry-breaking mass for $ |2| or by logarithmic 
radiative corrections . Our conclusions about the fine-tuning of the initial configuration that 
will lead to the onset of inflation are independent of the term providing this slope. The figures 
1-8 were generated for zero slope, and the constraints that we derived are strengthened if a 
significant slope is present. Including the imaginary parts of the scalar fields (some of which 
can be eliminated through appropriate gauge and i?-transformations along the D-flat directions) 
does not alter the qualitative picture. 

It must be pointed out, however, that our conclusions do not automatically apply to models 
with superpotentials different from that of eq. (p!8|). Models such as the smooth hybrid inflation 
of ref. 1^ may display less sensitivity on the initial conditions and their separate investigation 
is necessary. 

Acknowledgements: We would like to thank J. Garcia-Bellido, G. Lazarides and C. Panagio- 
takopoulos for many useful discussions. 
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Fig. 1: M = 10-^Mpi, A = 1, 5 = 1. 
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Fig. 2: M = 3 X lO^'^Mpi, X = l,g = l. 
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Fig. 3: Magnification of fig. 2. 
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Fig. 4: Magnification of fig. 2. 
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Fig. 5: M = lO-^Mpi, A = 1, 5 = 1. 
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Fig. 6: Magnification of fig. 5. 



0.8 



C\2 
O 



O 

I 

xo 

CO 

O 



0.6 



0.2 







4 <^ 



1 , 




1 

A 


o 


o 


i 

o 


A 


o 


1 

A 


o 


A 


1 

A 


o 


o 


1 

A 


A 


o 




A 


A 


o 


o 


o 


A 


o 


A 


o 


A 


A 


o 


o 


A 


A 


o 






o 


o 


o 


o 


A 


o 


A 


o 


A 


A 


o 


A 


A 


A 


o - 




O 


o 


o 


o 


o 


A 


o 


A 


o 


A 


O 


o 


A 


A 


A 


o 




O 


o 


o 


A 


o 


A 


A 


A 


o 


A 


O 


o 


A 


A 


A 


o 




o 


o 


A 


A 


o 


A 


A 


A 


o 


A 


O 


o 


A 


A 


A 


o— 




o 


o 


A 


A 


o 


O 


A 


O 


o 


A 


O 


o 


A 


A 


A 


o 




o 


o 


A 


A 


o 


O 


A 


O 


o 


A 


O 


o 


A 


A 


A 


o 




o 


o 


A 


A 


o 


O 


A 


O 


o 


A 


O 


o 


A 


A 


O 


o - 




o 


o 


A 


A 


o 


O 


A 


o 


o 


A 


o 


o 


A 


A 


o 


o 




o 


o 


A 


A 


o 


O 


A 


o 


o 


A 


o 


o 


A 


A 


o 


o 




o 


o 


A 


A 


o 


o 


A 


o 


o 


A 


o 


o 


A 


A 


o 








o 


A 


A 


A 


o 


A 


o 


A 


A 


o 


o 


A 


A 


o 


o 




A 


o 


A 


A 


A 


o 


A 


o 


A 


A 


o 


o 


A 


A 


o 


o 




A 


A 


A 


A 


A 


o 


A 


o 


A 


A 


o 


o 


A 


A 


o 


o - 




A 


A 


A 


A 


A 


o 


A 


o 


A 


A 


o 


o 


A 


A 


o 


o 




A 


A 


A 


A 


A 


o 


A 


o 


A 


A 


o 


o 


A 


A 


o 


o 




A 


A 


A 


A 


A 


o 


A 


o 


A 


A 


o 


o 


A 


A 


o 






A 


A 


A 


A 


A 


o 


A 


o 


A 


A 


o 


A 


A 


A 


o 


o 




A 


A 


A 


A 


A 


o 


A 


o 


A 


A 


o 


A 


A 


A 


o 


o 




A 


A 


O 


A 


A 


o 


A 


o 


A 


A 


o 


A 


A 


A 


o 


o - 




A 


A 


O 


O 


A 


o 


A 


o 


A 


O 


o 


A 


A 


A 


o 


o 




O 


A 


o 


O 


A 


o 


A 


o 


A 


O 


o 


A 


A 


A 


o 


o 




o 


A 


o 


o 


A 


o 


A 


o 


A 


O 


o 


A 


A 


O 


o 


o— 




o 


A 


o 


o 


A 


o 


A 


o 


A 


o 


o 


A 


A 


O 


o 


o 




o 


A 


o 


o 


A 


o 


A 


o 


A 


o 


o 


A 


A 


O 


o 


o 




o 


O 


o 


o 


A 


o 


A 


o 


A 


o 


o 


A 


A 


o 


o 


o - 




o 


O 


o 


o 


A 


o 


A 


o 


A 


o 


o 


A 


A 


o 


o 


o 




o 


o 


o 


o 


A 
1 


A 


O 


o 


A 


o 


o 


A 


A 


o 


o 


o 


1^ 


(-1 








ft 






A 




J, 


A 


n 













0.2 



0.4 0.6 

lO^xa 



0.8 



Fig. 7: Magnification of fig. 5. 
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Fig. 8: M = 10~'^Mpi, A = 1, 5 = 1. Initial total energy density E = 10~"^M^; and (j) = a > 



